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The space of positive definite symmetric matrices has been stud- 
'/^ , ied extensively as a means of understanding dependence in multi- 

variate data along with the accompanying problems in statistical 
inference. Many books and papers have been written on this sub- 
L^ , ject, and more recently there has been considerable interest in high- 
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dimensional random matrices with particular emphasis on the distri- 
bution of certain eigenvalues. With the availability of modern data 
1 -^ . acquisition capabilities, smoothing or nonparametric techniques are 

j^ ' required that go beyond those applicable only to data arising in Eu- 

clidean spaces. Accordingly, we present a Fourier method of minimax 
Wishart mixture density estimation on the space of positive definite 
symmetric matrices. 

—^ ' 
> 
^Sj . 1. Introduction. The space of positive definite symmetric matrices has 

■^ \ been studied extensively in statistics as a means of understanding depen- 

dence in multivariate data along with the accompanying problems in sta- 
tistical inference. Many books and papers, for example, [7-9, 17, 19, 22] 
and [23], have been written on this subject, and there has been consider- 

^^ . able interest recently in high-dimensional random matrices with particular 

emphasis on the distribution of certain eigenvalues [11] and on graphical 
models [15]. 
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In this paper we consider the problem of estimating the mixing density of 
a continuous mixture of Wishart distributions. We construct a nonparamet- 
ric estimator of that density and obtain minimax rates of convergence for 
the estimator. Throughout this work, we adopt, as a guide, results developed 
for the classical problem of deconvolution density estimation on Euclidean 
spaces; see, for example, [2, 4, 5, 14, 18] and [26]. Much of the difficulty 
with the space of positive definite symmetric matrices is due to the fact that 
mathematical analysis on the space is technically demanding. Helgason [10] 
and Terras [25] provide much insight and technical innovation, however, and 
we make extensive use of these methods. 

We summarize the paper as follows. In Section 2 we discuss and set up 
the notation for Wishart mixtures. In Section 3 we begin by reviewing the 
necessary Fourier methods which allow us to construct a nonparametric 
estimator of the mixing density, and then we provide the estimator. The 
minimax property of our nonparametric estimator is stated in Section 4 
along with supporting results. Section 5 presents simulation results as well 
as an application to finance examining real financial data. Finally, Sections 6 
and 7 present the proofs. 

2. Wishart mixtures. Throughout the paper, for any square matrix y, we 
denote the trace and determinant of y by try and |y|, respectively; further, 
we denote by I^ the mxm identity matrix. We will denote by Vm the space 
of m X m positive definite symmetric matrices. 

For s = (si, . . . , Sm) G C*" with Re{sj + • • • + Sm) > (j - l)/2, j = 1, . . . , m, 
the multivariate gamma function is defined as 

m , 

(2.1) r^(si,...,s^)=7r'"('"-^)/^nr{^. + --- + «--2(j'-l) 

where r(-) denotes the classical gamma function. 

We denote by G the general linear group GL(77i,M) of all m x m nonsin- 
gular real matrices, by K the group 0(?tt,) oi rax m orthogonal matrices 
and by A the group of diagonal positive definite matrices. The group G acts 
transitively on Vm by the action 

(2.2) GxV^^Vm. {g,y)^g'yg: 

g € G, y € Vm, where g' denotes the transpose of g. Under this group 
action, the isotropy group of the identity in G is K; hence the homoge- 
neous space K \G can be identified with V^n by the natural mapping from 
K \G ^ Vm that sends Kg i— >■ g' g. In distinguishing between left and right 
cosets, we place the quotient operation on the left and right of the group, 
respectively. 
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For y = [yij) G Vm, define the measure 

It is well known that the measure d-^y is invariant under action (2.2). Rela- 
tive to the dominating measure d*y, the probability density function of the 
standard Wishart distribution with A^ degrees of freedom is 

('•') "(^^ = 2^W2r^(o|...,0,Ar/2) '^'''^'^"p("^''' 

y G Vm- Consequently, for a G Vm, we note that tr(cr~^' ^yo""^' ^) = tT{a~^y) 
and \cj~^''^y(T~^''^\ = \a~^y\. It then follows that, relative to the dominating 
measure d^,y, the density of the general Wishart distribution, with covariance 
parameter a, is w{a~^y), y G Vm- 

Suppose next that u is a random matrix and, relative to the dominating 
measure d*(T, has a continuous mixing density, /, that is invariant under 
the action (2.2). By integration with respect to a, the continuous Wishart 
mixture density is given by 

(2.4) r(y)= / f{a)w{a~'y)d,a, 

JVra 

y G Vm- For the case in which m = 1, the standard Wishart density is es- 
sentially a chi-square density, in which case (2.4) is a continuous mixture of 
chi-square densities. 

In general, (2.4) is a convolution operation for functions on Vm- We de- 
note by x^''^ any matrix with x''''^x^''^ = x, where x*' ^ = {x^'"^)' and denote 
x"*'^ = (x~^'^)*. Define X o Z, the convolution of two random matrices X 
and Z which are distributed on Vm by 

and /i * /2, the convolution of /i G L^{Vm) and /2 G L^iVm) by 

(/l * /2)(y) = / /l(x)/2(x-*/2yx-l/2)d,X for y G Vm, 

where L'i{Vm) is the space of integrable functions raised to the g'th power 
on Vm for g > 1. If X and Z with densities / and w, respectively, are inde- 
pendent, then Y = X o Z has the density r = f *w since w{a~*''^ya~^''^) = 
w{a~^y). Finally, (2.4) can be transformed into a scalar multiplication; see 
Section 3.2. 

3. Fourier analysis on T^^n and estimation of the mixing density. In this 
section we review the Fourier methods needed to transform the convolution 
product (2.4) and to construct a nonparametric estimator of the mixing 
density /. 
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3.1. The Helgason-Fourier transform. For y € Vrm denote by \yj\ the 
principal minor of order j, j = 1,. . . ,m. For s G C"^, the power function ps : 
V.m ^ C is 

m 

(3.1) Ps{y) = ll\yj\'\ 

y € Vm- Let d^/c denote the Haar measure on K, normahzed to have total 
volume equal to one; then 

(3.2) h,{y)= I ps{k'yk)d,k, 

Jk 

y € Vm, is the zonal spherical function on Vm- It is well known that the 
functions hs are fundamental to harmonic analysis on symmetric spaces 
[10, 25]. If si, . . . ,Sm are nonnegative integers then, up to a constant factor, 
(3.2) is an integral formula for the zonal polynomials which arise in many 
aspects of multivariate statistical analysis [19], pages 231 and 232. 

Let C^{Vm) denote the space of infinitely differentiable, compactly sup- 
ported, complex- valued functions / on Vm', also, let 

C^iVm/K) = {/ G C^iVm) : fik'yk) = f{y) for ah A: G K, y € Vm}- 

For s G C" and k G K, the Helgason-Fourier transform ([25], page 87) of 
a function / G C^iVm) is 



(3.3) nf{s,k)= f{y)psik'yk)d,y, 



where ps{k'yk) denotes complex conjugation. 

For the case in which / G C^{Vm/K), we make the change of vari- 
ables y I— 7> k[yki in (3.3), ki G K, and integrate with respect to the Haar 
measure d^/ci. Applying the invariance of / and formula (3.2), we deduce 
that 7if{s,k) does not depend on k. Specifically, 7if{s,k) = f{s) where 



(3.4) f{s)= f{y)hsiy)d.y, 

JVra 

s G C™, is the zonal spherical transform of /. 

In the case of the standard Wishart density (2.3), which is a i('-invariant 
function, the zonal spherical transform is well known (Muirhead [19], pa- 
ge 248; Terras [25], pages 85 and 86): 

..X rm(Sm-l,...,Sl,-(siH \-Sm)+N/2) (l 

wis) = — — r^ n,o — Ir 

^ ' r^(0,...,0,A^/2) ^V2 

3.2. The convolution property of the Helgason-Fourier transform. The 
following result shows that the convolution operation can be transformed 
into a scalar multiplication. 
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Proposition 3.1. Suppose X and Z with densities fx and fz, re- 
spectively, are independent, and Z is K -invariant. Let /y be the density 
ofY = XoZ. Then 

nfY{s,k)=nfx{s,k)fz{s) forseC"' andkeK. 

Proof. Note 

■HfY{s,k) =E,ps{k'Yk) =Eps{k'X^/'^ZX'^/'^k). 

Using the iT^A'^-Iwasawa decomposition of X^'^k (Terras [25], page 20), we 
have X^''^k = HU for H £ K and U, an upper triangular matrix. Observe 

Ep-s{k'X'/^ZX^/^k) = Ex{p-siU'U)EzixPs{H'ZH)} 
= fz{s)Ex{PsiU'U)} 
= nfx{s,k)fz{s), 
where Proposition 1 of Terras [25], page 39, is used for the first equahty. D 

3.3. The inversion formula for the Helgason-Fourier transform. For oi, 
a2 € C with Re(ai),Re(a2) > 0, let 

r(ai)r(Q2) 

B[ai,a2) = —( — - — T 
r(ai + a2) 

denote the classical beta function. For s € C"* such that Re(sj + ■ ■ ■ + Sj)> 
— gO — i + 1) for all l<i < j <m — \, the Harish-Chandra c-function is 

B{l/2,s^ + --- + Sj + {j-i + l)/2) 

l<i<j<m—l 

Let p= (|, . . . , ^, j{l — m)), and set 



^^^ ""^^'.JL . Bil/2,{j-i + l)/2) 



nr=ir(j72) 



(3-6) '^™"(27ri)-W-+i)/4m!' 

(3.7) C™(p) = {seC"^:Re(s) = -/j} 

and 

d*S = Ujn\Cm{s)\~'^ d*Si ■ • ■ d*Sm- 

Let M = {diag(ibl, . . . , ±1)} be the set of m x m diagonal matrices with 
entries ±1 on the diagonal; then M is a subgroup of K and is of order 2™. 
By factorizing the Haar measure d^/c on K, it may be shown ([25], page 88) 
that there exists an invariant measure d^k on the coset space K/M such 
that 

f dJ = L 

Jk&K/M 
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The inversion formula for the Helgason-Fourier transform % in (3.4) is 
that if / G C^iVm), then [10, 25] 



(3.8) f{y)= / 'Hf{s,k)ps{k'yk)d,~kd.s, 

y G Vm- In particular, if / G C^{Vm/K), then 



fiy)= / fis)hsiy)d^s, 

y G Vrm and there also holds the Plancherel formula, 
(3.9) / |/(y)|2d,y= / / \'Hf{s,k)\^dJd.s. 

JPm JC'"(p) JK/M 

We refer to Terras [25], page 87 ff., for full details of the inversion formula 
and for references to the literature. 

3.4. Eigenvalues, the Laplacian and Sobolev spaces. For y = {yij) G Vm, 
we define the m x m matrix of partial derivatives, 

.^ + 0i 



dy \2 dyij 

where 5ij denotes Kronecker's delta. The Laplacian, A, on Vm can be written 
([25], page 106) in terms of the local coordinates yij as 

A = -tr( (y 



dy, 

The power function ps in (3.1) is an eigenfunction of A (see [19], page 

229, [21], page 283, [25], page 49). Indeed, let rj = Sj + Sj+i H h s^ + 

^{m — 2j + 1), j = 1, . . . , m, and define 

(3.10) A, = -{rl + ... + rl) + ^mim^ - 1); 

then Aps(Y) = XsPs{Y). Since Re(s) = —p then each rj, j = 1, . . . ,m, is 
purely imaginary; hence, A^ > 0, s £ Cm{p)- 

The operator 71 changes the effect of invariant differential operators on 
functions to pointwise multiplication: if / G C^{Vrn), then 

n{Af){s,k) = Xsnf{s,k), 

s G C™, k £ K ([25], page 88). For c/5 > 0, we therefore define the fractional 
power, A'^'^, of A, as the operator such that 

H{A'^'^f){s,k) = \tl^Hf{s,k), 

f G C^iVrn)- Having constructed A*^'^, we define the Sobolev class, 

-F^ = {/gC-(P„):||A^/2/||'<oo}, 
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where for f£C°°(Vm), 

1/11 = (/ \f{y)\'d.y' 

denotes the L^('Pm)-norm with respect to the measure d,,y. For Q > 0, we 
also define the hounded Sobolev class, 

J-^(Q) = {/ G C^{V^) : IIA'^/^f < Q}. 

4. Main result. In this section we wih present the main result. We do so 
by applying the Helgason-Fourier transform to the mixture density (2.4) so 
that 

(4.1) nr{s,k) = nf{s,k)w{s), 

s € C", k £ K; see Proposition 3.1. Having observed a random sample Yi, 
. . . ,Yn from the mixture density, r, in (2.4), we estimate 1-Lr{s,k) by its 
empirical Helgason-Fourier transform, 

I "^ 

(4.2) nnT{s,k) = -Y^Vs{k'Y,k). 

£=1 

On substituting (4.2) in (4.1), together with the assumption that w{s) ^ 0, 
s S C™', we obtain 

n.fis,k) = ^h^, 

w{s) 

sGC™, A;GK. 

Analogous with classical Euclidean deconvolution, we introduce a smooth- 
ing parameter T = T{n) where T{n) — )■ oo as n — > cxj, and then we apply the 
inversion formula (3.8) using a spectral cut-off based on the eigenvalues of A. 
First, we introduce the notation 

C"\p,T) = {seC"\p):Xs<T}, 

where C'"(p) is defined in (3.7). We now define 

(4.3) Uy)=f I ^''lf:^\ s{k'yk)dikd,s, 

y G Vm, and take this as our nonparametric estimator of /. 

We now state the minimax result for the estimator (4.3). Let C denote 
a generic positive constant independent of n. For two sequences of real num- 
bers {an} and {bn}, we use the notations a„ <^ bn and a„ » 6„ to mean 
On < Cbn and On > Cbn, respectively, as jz — >■ oo. Moreover, a^ x bn means 
that On <C bn and a„ ^ bn ■ 
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Theorem 4.1. Suppose f is a density on Vm and N > {m — l)/2. Then, 
for the Wishart mixture (2.4), 

(4.4) sup E||/,-/f «(logn)-2v' 

and for any estimator gn of f , 

(4.5) inf sup E||5„-/||2»(logn)-2'^. 

We now provide some comments about this result. In the situation whe- 
re (2.4) is a finite sum, so that 

q q 

r{y) = '^fiw{aj^y) and ^/£ = 1, 
i=i e=i 

we have the finite mixture model. Methods for recovering the mixing coeffi- 
cients can be covered by the techniques employed in [3]. We note that the 
continuous mixture model is a generalization of this approach. 

It is noted that the condition / is a density and seems to be mild. The 
upper bound property of (4.4) is established in [13], Theorem 3.3, with 
/3 = 1/2. In the latter, the moment condition 

(4.6) / |yi|-i..-|y^_iri|y|("-^)/V(y)d.y<oo, 

on the principal minors |yi |, . . . , \ym\ of y G Vm is assumed. In our theorem, 
we did not impose this moment condition as condition (4.6) is automatically 
satisfied. This is pointed out and commented upon below in the proof. 

To derive the lower bound for estimating / in the L^('Pm)-norm, we shall 
follow the standard Euclidean approach. Thus we choose a pair of func- 
tions f^, /"" € T^{Q), and, with w denoting the Wishart density (2.3), we 
shall show that, for some constants Ci,C2 > 0, 

Iir-/°f >Ci(logn)-2'^ 

and 

(4.7) ^2^f^u^jrr^u^)<9l^ 

n 
where 

X [91,92)= / r-. d*y. 

Precisely, let us suppose we can choose f^ S J^,p{Q) and a perturbation 
ip € Fip{Q), and, for 5 = 5n> ^, let ip be a scaling of if) such that ||V' || >< 
5-^/'^U\\- Define 



/" = /o + Cv,r^+^/V- 
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If 6 can be chosen so that 

x\f*w,r*w)<Cn-\ 

then the lower bound rate of convergence is determined by f^"^*^. We shall 
develop such a construction and, moreover, do so in a way such that 6 x log n 
as n — )• oo. 

Remark 4.2. The profound influence of Charles Stein on covariance 
estimation originates largely from his Rietz lecture; see [22]. The idea is 
that for certain loss functions over Vm , the usual estimator of the covariance 
matrix parameter is inadmissible. Through an unbiased estimation of the 
risk function over covariance matrices. Stein was able to improve upon the 
usual estimator by pooling the observed eigenvalues of the sample covariance 
matrix. Subsequent to this, through a series of papers, improvements were 
obtained in Half [7-9]. Other related works include Takemura [24], Lin and 
Perlman [16] and Loh [17], to name a few. 

In this paper, we contribute to the case in which one observes data from 
a continuous mixture of Wishart distributions, not merely a sample from 
a single distribution. Therefore, the parameter of interest would be the mix- 
ing density of the covariance parameters. And the nonparametric estimator 
of the mixing density (4.3) is an attractive candidate because of its minimax 
property. Based on this procedure, one could consider the moment, or mode, 
of fn, as a possible estimator of the corresponding population parameters. 
Alternatively, one could take a nonparametric empirical Bayes approach as 
inPensky [20]. 

5. Numerical aspects and an application to finance. This section presents 
numerical aspects for the m = 2 case with an application to finance. 

5.1. Computation of estimators. Suppose X and Z are independent 
with Z having a Wishart distribution. Let Y = X^I'^ZX^I'^ where X^/^ is 
upper triangular. For visualization, we display estimators of the marginal 
density for D where X = H'DH with H eK and D e A. Let 



fn{s) = -yZhs{E,) and fnis'-''''^'^ 

71 ^ ^ 



n 



w{s) ' 



where Ej G A_|_ denotes the diagonal matrix of eigenvalues of Y^-, j = 1, . . . , n. 
Denote by /^ the density of eigenvalues of X. Then, the estimator for f^ 
is given by 

(5.1) fn{a)= Re{fn{s)hs{a)}d:,s for a = diag(ai,a2) G ^. 

Consider the computation of hs{a) when 

s = -p + ib = (-1/2 + ibi, 1/4 + ib2) 
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SO that Re(s) = —p. From pages 90 and 91 of [25], the spherical function is 
given by 

where Legendre function P-i/2+it{x) can be computed using conicalP_0 (t , 
x) in the gsl hbrary in R. 

The mean integrated squared error (MISE) of /„ is defined by 

MISE(/„) = E / (/f (a) - /^(a))' d,a. 
Ja 

It is reasonable to choose T which minimizes MISE(/^) or equivalently 



M{T)=E / {f^{a)yd,a-2E / /„^(a)/^(x)d,a. 
Ja Ja 

One can find an unbiased estimator Mq(T) of M{T); see [12]. We choose T 
by 

arg min Mq (T) . 

T 
Monte Carlo approximation is used for integration of (5.1) and Mq{T). 

5.2. Simulation. Denote by WN{a) the Wishart density with degrees 
of freedom N and covariance matrix a. We generate data as follows. For 
j = l,...,n: 

• generate Zj ~ W2o(/2); 

• generate Xj i^ f; 

• do a Cholesky decomposition of Xj = (Xj)*'^(Xj)^/^, and calculate Yj = 

{x,r/^z,{x,f/\ 

As examples, we consider a unimodal mixing density ^^15(2/2) and a bi- 
modal density 0.5Vl^i5(2/2) + 0.5^15(6/2). Figure 1 show the results for the 
unimodal case whereas Figure 2 show the results for the bimodal case. In 
each of these plots the domain consists of the two eigenvalues starting with 
the largest. One can see that the general shapes of the estimators become 
closer to that of the true density as n increases. 

5.3. Application to stochastic volatility. Stochastic volatility using the 
Wishart distribution is of much interest in finance; see, for example, [1] 
and [6]. In particular, this entails a situation precisely of the form (2.4). 
Let us apply this to the situation where we are interested in estimating the 
mixing density. 

Although our methods can be applied to a portfolio of many assets, let us 
restrict ourselves to two assets since this would be the smallest multivariate 
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Fig. 1. Unimodal case: upper left displays the true density 0/^15(2/2), upper right shows 
an estimate with n = 500, lower left with n = 1,000 and lower right with n = 2,000. 



example. Indeed, let 5j and 5| denote the daily closing stock prices of 
Samsung Electronics (005930. KS) and LG-display (034220.KS), respectively, 
traded on the Korea Stock Exchange (KSC) for 2010, where the data can 
be easily accessed on public financial websites. We will assume as usual that 
Q^ = log{S^_^i/S^) follows a bi-variate normal distribution for A; = 1,2. We 
transform the daily data to weekly data and compute the weekly 2x2 
covariance matrix Yi for i = 1,...,52. In case a week has a hohday, we 
repeat the last previous observation. Under the usual assumptions this would 
constitute observations from a mixture model (2.4) with a standard Wishart 
distribution with four degrees of freedom. 

Figure 3 plots the mixing density estimator corresponding to the two 
eigenvalues. One can see that there are two peaks, suggesting a possible 
bimodal stochastic volatility mixing density in the eigenvalues. 



12 
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Fig. 2. Bimodal case: upper left displays the true density 0.5^15(2/2) + 0.5^15(6/2), 
upper right shows an estimate with n = 500, lower left with n = 1,000 and lower right with 
n = 2,000. 



6. Proof of upper bound. The strategy here is, first, to decompose the 
integrated mean-squared error into its variance and bias components, 

E|l/„ - ff = E||(/„ - E/„) + (E/„ - /)f 
(6.1) 

= E||/„-E/„f + ||E/„-/||2, 

and, last, to estimate each component separately using estimates based on 
the Plancherel formula and the inversion formula for the Helgason-Fourier 
transform. 

6.1. The integrated bias. 

Lemma 6.1. Suppose that f G J^^{Q) and (p > dim'Pm/2. Then 

mu-ff^T'^. 
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Fig. 3. A density estimator for the weekly covariance matrix of stock prices. 



Proof. We have for x G P, 

Efn{x) - fix) 



C^{p,T) Jk<^K/M 



T-if{s, k)ps(k'xk) d*/c d*s 



(6.2) 



C^Cp) Jk&K/M 



V.f{s, k)ps{k'xk) d^fk d^=s 



' As>T,Re(s)=-p Jk£K/M 
Applying the Plancherel formula, we obtain 

mu-ff 



T-Lfis, k)ps{k'xk) d*A;d*s. 



Xs>T,Kc(s)=~pJK/M 



\Hf{s,k)\^d^kA^s. 



Consequently, 

l|E/n-/||' 

(6.3) 



\s>T,Kc{s)=-p JK/M 



1-^/(8, ^)|2d,fcd*S 



'Xs>T,Rc{s)=-p Jk/M 



<T~'P / A^|?^/(s,fc)|2djd*s, 

Jc'^ip) Jk/m 

where we use the fact that 

xt\nf{s,k)\^^\xf''nfis,k)\^ = \'H{A'^/^f){s,k)\^ 
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Therefore 



J<C^(p) JkeK/M 



where the equahty follows from the Plancherel formula. By assumption, 
/ S T^{Q), the latter integral is bounded above by Q, so we obtain 

mfn-ff<QT-^, 

and the proof is complete. D 

6.2. The integrated variance. To obtain bounds for the integrated vari- 
ance, several preliminary calculations are needed. In particular, we begin 
with the variance calculation of the empirical Helgason-Fourier transform, 
which has similarities to the usual empirical characteristic function. 

Lemma 6.2. For s e C^ip) and k e K/M, 

^■Hnr{s,k) -^■Hnr{s,k)\^ = -{\Hr{-2p,k)\^ - \-Hr{s,k)\^). 

Proof. By (4.2), 

\'Hnr{s, fc)p = Hnris, k)'Hnr{s, k) 

(6.4) =^y\ Ps{k'Y,k)ps{kYek) 

= ^1 E \Psik'Y,k)\^ + Y.Ps{k%k)ps{k%k) I. 






Observe also that 



Ps{w)ps{w) = \wi\''^^ ■ ■ ■ |i^m|*'"|^'^l|^^ ■ ■ -{W: 



m\ 



- U,,j2Re(si).. I |2Re(s„) 

— I "'I I I "Jm\ 

= p-2p{w) 

since Re(s) = —p. Applying this result to (6.4) and taking expectations, we 
obtain 

E\nnris,k)\^ = ^eI Y, \Psik'Y,k)\^ + Y,Ps{k'Yjk)p,{k%k) I 
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= ^1 X^Ep_2p(A;'y,A:) + J^Ep,(fc'y,fc)Ep,(A:'y,A;) i 

= -nr{-2p, k) + !^|^r(s, fc)|2, 
n n 

where the last equahty follows from the fact that Yi, . . . ,Yn are independent 

and identically distributed as Y, and because Kps{k'Yk) = 'Hr{s,k). D 

Following Terras [25], pages 34 and 35, let 

A+ = {a = diag(ai,...,am) G^:ai > •■• > flm} 
denote the positive Weyl chamber in A. For a = diag(ai, . . . ,am) S A, let 
da = YYjLi ^7 ^"^i s-'^d set 

m 

and define the normalizing constant bm by 6~^ = vr"'"* ~^"^'''^Y[T=iJ'^U/'^)- 
Denote d^,a = 6^7(0) da. 

Lemma 6.3. ^sr^>oo, 

jidimPm/S 

E||/„-E/„||2< sup |u'z(s)r' . 

seC'"(p) '^ 

Proof. By the Plancherel formula, 

E||/„-E/„f 



= / / E\'Hnfis,k)-EHnf{s,k)\^d~kd,s 

JC™(p,T) Jk/M 

= [ I E\'Hnr{s,k)-Ennr{s,k)\^Ak\wz{s)\"'^d^s 

Jc"^{p,T) Jk/M 

<- sup |u;z(s)rM \nr{-2p,k)\^ dk d*s 

^ \.,<T,Rc(s)=-p Jk/M Jc"^(p,T) 



^ As<t,Rg(s)=-p Jk/M jc™(p,t) 

2idim'Pm/2 

^ sup \wz{s)\~'^ 

Xs<T,Rc{s)=-p "■ 

as r — >■ cx) . 

Choose a € A~^ . Observe that 
(6.5) p.,,{a) = «-(™-i)+(™"i)/2 . . . ^-i+(™-i)/2^(,n-i)/2 < 1. 

Since p-2p{k'ak) is a continuous function of /c on a compact set K, p^2p{k'ak) 
is uniformly bounded on K such |p_2p(/c'a/c)| < C on i^. Since / is a density 
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SO that r is also a density, we have 

\'Hr{-2p,Im)\= / r(y)p_2p(y)d*y 

■JVm 



r{k' ak)p^2pik' ak) d*ad*A; 
A+ JK 



< I I r{k'ak)\p-2p{k'ak)\d^ad^k 
Ja+ Jk 

<C / / r{k'ak)d^ad^k 
Ja+ Jk 



= C r(y)d*y 

JVm 

= c. 

Hence, it follows from continuity and the compactness of K/M 



I \Hr{-2p,k)\^dk<oo, 
Jk/m 



I K/M 

which has been used in the above calculation. 
In addition, we use the fact that, as T — )• cx), 

sup |c^(s)|-2«r"^(— 1)/^ 

C™(p,T) 

a result which follows from Proposition 7.2 of Helgason [10], page 450. D 

The proof of the upper bound can now be obtained by applying Lem- 
mas 6.1 and 6.3 to (6.1) and setting Tx (logn)^. 

7. Proof of lower bound. We need to provide some detailed calculations, 
and the essence of the proof is contained for the case m = 2; hence we will 
keep this assumption for the remainder of this paper. The generalization to 
m> 2 may be obtained by using higher order hyperbolic spherical coordi- 
nates. In this section, we assume that ■!/; is a i('-invariant function defined 
on V2- 

7.1. Convolution and Helag son-Fourier transform in polar coordinate. 
For y e 7^2) let y = k'ak with a = diag(ai, 02) G A'^, k ^ K so that ipiy) = 
ip{a). Let 

a = Du,e''^ 

with Dz = diag(e^, e~^) for z € M, and write 
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By a change of variables, 

(7.1) / ip{y)d^y= / VMd*n, 

where 



and 



P = {n : til e M+ and U2 G M} 



d*ii = 47rsinhui dui du2. 



Denote kg = [gj^g cos6» ]■ '^^^ next lemma is straightforward so we shall 
omit the proof. 

Lemma 7.1. For ui G M"*", vi G M"*" anc? 6 G [0, 27r], i/ie matrix equation 

k^Dnk'^ = D^^^/2keDv,k'QD^^^/2 

has a solution R* = R*{ui,vi, 9) and ^* = ^*{ui,vi,9). Further, coshi?* has 
minimum and maximum values cosh(ui — vi) and cosh(tii + t;i), respectively, 
and R* and ^* can be defined uniquely. 

In general, if both / and g are ET-invariant functions on Vm, then f * g is 
also iT-invariant and f * g = g * f- Hence, ip *w is i('-invariant and ip *w = 
w * ip due to i^-invariance of ip and w. From this and Lemma 7.1, we can 
define 

^u,vA^) = ^ f\{R*{e,ui,vi),z)de 

for z G M, til G M"^, vi G K^. Denoting by W the distribution function cor- 
responding to the standard Wishart density w{u) with respect to the mea- 
sure d^,ti, we have that for v (^ D, 

(7.2) {^p^nj)iv)= [ ^u,vAv2-U2)dWiu). 

Jv 

The Laplacian for iC-invariant functions in polar coordinate is given by 

A = A^^ -I- Au2 , where 

d 92 92 

A„,=-coth(tzi)— -^, A„, = -^, 

and the spherical function is given by 

(7.3) /is(n) = P,^ (cosh til )e(^i+2s2)«2 for n G P 

with Ps, the Legendre function; see Terras [25]. It can be seen that 

Ani/ls = -ls{s + l)hs, Auihs = -f (si + 2s2fhs, 

so that A/is = Xghg with As = — ^{(si + 2s2)^ + s{s + 1)}. 
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The Helgason-Fourier transform of ip is given by 



(7.4) V'(s)= / 'ip{u)Ps,icoshui)e^'^+^'^'>''^d^u 

Jv 

from (7.1) and (7.3). Suppose ip is separable so that ip{u) = "01(^1 )V'2(^t2) ■ 
Then, (7.4) imphes 

(7.5) i;is)=MMsi)^Msi + '^S2), 

where C and M denote, respectively, the Laplace transform and the Mehler- 
Fock transform; see Terras [25]. 

7.2. x^ -divergence. Choose ipi as the perturbation function as in Fan [5]. 
Then, one can construct V^(u) satisfying the following conditions: 

(PI) ip is ET-invariant and separable with ip{u) = ^'1(^1)^2(^2) for u ^ D. 

(P2) Cij2iit) = 0fort^[l,2]. 

(P3) ^gJ-^(Q). 

(P4) Vi(ui) = O(cosh-"'0'Ui) and 02('U2) = O(e""'°l"^l), where 0<C<1 

and mo^>(iV- 1) (1-0/2- 
(P5) f^^{u)d^u = 0. 

Let ph be a density on M such that ph is sufficiently smooth and satisfies 
Pb{u2) = Cfeexp(— 6|u2|), \u2\ > Co, where c^ is a normalizing constant. Define 
the function 

/0(n) = C6(coshni)~Vb('U2), 

where Cb = Cb{b - l)/(27r) for 6 > 1. 

For a function §-.1^2^^ and 5 > 0, define 

/(y) = 5(|y|^'"'^/'y) foryeV2 

so that 

ipiu) =ip{ui,5u2) foruGV. 



Proposition 7.2. Suppose (P1)-(P5) hold. For a pair of densities 

f and /" = /0 + Cv,r'^+^/V, 
the x^ -divergence between g^ = f^ * w and g"^ = f^ * w satisfies 

x\9',9n<C/n 
provided that b < 2 niin(37r, (A'' — 1)(1 — ^ — 1). 

The fact that (logn)"^*^ is a lower bound follows from Proposition 7.2 
whose proof follows from a sequence of lemmas below. 
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7.3. Perturbing function. Denote ((71,92) = (si,si + 2s2) and /3j = Im(gj) 
for j = 1, 2. Note that ^2 = «/32 for s e C2(/3). 

Lemma 7.3. Suppose (PI) /lo/ds. T/ien ||-0'^|| = (5~^|1^||. 

Proof. By (7.5) and the change of variable U2 i— )• Su2, we obtain 

The desired result follows from the Plancherel formula (3.9) and change of 
variable s >-^ q. D 

Lemma 7.4. Suppose (PI), (P2) and (P3) hold. Then there exists a posi- 
tive constant C^ such that C^(5~'^+^'^'0^ € Fip{Q). 

Proof. Suppose s G C^(/)) n S. Observe that 

A. = -\{si{si + 1) + (si + 2s2f} = \{Pl + /3| + \) 
and that for (5 > 1, 

Now, the Plancherel formula (3.9) and change of variable s^^ q gives 

A suitable choice of C^, gives the desired result. D 

Lemma 7.5. Under (PI) and (P2), HV'^ * wjf < C^e-^'^^HV^f. 

Proof. For s G C2(p) n 5, 

|u)(s)P < Ce-''("i+2"2) ^ c'e-2'r/32 < Ce''^''. 

The desired result follows from the inequality ^-'^i^^+'^'^V" < Ce~^'^ , the 
Plancherel formula (3.9) and change of variable s i— )■ g. D 

7.4. Tai/ behavior. 

Lemma 7.6. For ci > and C2,C3 gM, 
/ dW{u)=A-K exp{(iV-l)ti2-e"2coshci}dn2. 

J {ui>Ci,C2<U2<Ca} 'J{C2<U2<C3} 

Proof. Change of variable gives the desired result. D 



20 HAFF, KIM, KOO AND RICHARDS 

If U has the distribution function W, then {Ui,6U2) has the distribution 
function Ws with density 6~^w{ui,U2/6). By (7.2) and a change of variables, 

(7.6) {iP^ * W){V) = [ ^u,vA^V2 - U2)dWs{u). 

Jv 



Lemma 7.7. We have g^ = f^ *w is K -invariant, and as vi ^- oo and 

1^2! — > C)0 

5°(^)>Ce-*(''i+l^2|)_ 



Proof. Choose a constant C1/2 such that 

dW{u) > 1/2 



'{0<ni<Ci/2,|n2|<Ci/2} 

andph{u2) = Cbexp(— 6|n2|) for JU2J >Ci/2- The desired result follows from (7.2) 
and Lemma 7.1. D 

Lemma 7.8. Define 

r]s{v) = I cosh-'""(ni - vi)e-'^°\''^-''^\ dWs{u). 
Jv 

Then, there exist M and a constant C such that when vi > M and \v2\/S > M , 

m{v) < Cexp{-i(iV - 1)(1 - 0(^1 + \v2\/6)} 
for all 6>l provided that < ^ < 1 and moC > {N - 1)(1 - ^)/2. 

Proof. Denote 

Jii = {u:\ui -Vi\ <£.'V1,\V2-SU2\ <^|f2|}, 
J12 = {u:\ui -Vl\ <C'V1,\V2-SU2\ >C\V2\}, 
J21 = {u:\ui -vi\ >Cvi,\v2-Su2\ <C\V2\}, 

J22 = {U- \U1 -Vl\ > ^Vl, \V2 - dU2\ > ^\V2\} 

and for i,j = 1,2, 

lij = I cosh-'"«(ni - yi)e-"^ol''2-"2l dWs{u). 

Consider In. Suppose ^2 < 0. If |ui — wi| < ^vi and \v2 — 6u2\ < ^1^21) then 

(1 - 0(^1 + \v2\/5) <Ui-U2< (1+0(^^1 + \v2\/5)- 

Let C = (1 - 0{vi + \v2\/5) and A = e^. Note that 



111 < C ( / dW{u) + / dW{u) ) := I^ + /2 
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Since e"^ cosh(ti2 + C) = ^(-^^^"^ + A^^), then we have 

J{0<t<l/X} 

and by Lemma 7.6 with ci =U2 + C, C2 = — C ^-iid C3 = 

l2=c[ e(^-^)"2 exp(-e"2 cosh(u2 + C)) dna = ©(A^^^-^^/^^. 

Hence, if ^2 <0, then In = 0{X~^^~^^''^) as fi + |u2|/(^— s-oo. Suppose t)2 > 0. 
If \v2 — du2\ < ^1^2! and \ui — vi\ < ^fi, then 

(1 - 0{vi + \V2\/S) <Ui+U2< (1 +6('^1 + \V2\/S). 

Observe that 

Iu<c( f dW{u) + / dWiuU := 1+ + 1+. 

\J{mi>0,«2>C} J{mi>C-M2,0<M2<C} / 

Since e"^/^ cosh(C - U2) = ^(A + A'^e^^^)^ ^hen we have 

/+ = 47r / t^-2e-* dt = 0{e^^) 
J{t>\} 

and by Lemma 7.6, with ci = U2 — Ci C2 = and C3 = C, 

4 = C I e^N-i)u2 exp(-e«2 cosh(C - ^2)) dua = ©(e"^). 

i{0<«2<C} 

Hence, for ^2 G M, 

In = 0(A-(^-i)/2) = 0(exp{-i(iV - 1)(1 - 0(i^i + \v2\/S)}) 

as fi + |f2|/(5 — >■ 00. 

Consider I21. Let C = (1 — Ol^al/^ and A = e^. Suppose ^2 < 0. On J21, 
-^2 < -(1 - 01^2!/'^ = -C on J21. Note that cosh(^ui)-™o < (ie^^'i)"""^ < 
2-moe-(i/2){iV-i)(i-0^i. If fohows from these and Lemma 7.6 that 

/21 = 0(exp{-i(iV - 1)(1 - e)C"i + \V2\/S)}), 

and if V2 > 0, 

hi = 0(exp{-i(iV - 1)(1 - e)(«i + |t'2|/<5)}). 

Observe that mQ^\v2\ > mo^\v2\/6 > ^{N — l){l — ^)\v2\/6. It fohows from 
this and Lemma 7.6 that 

)5b2| 

Iu2&RJui>{1-^)vi_ 



I12 < Ce-'"o«l"2| f I dVF(tx) 



0(exp<!-^(iV-l)(l-e)(t;i + b2|/5) 
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and 

/22 < Ce-"'o«l^'2l(cosh(et;i))-'"" / dW{u) 

Jv 

= o(expl-^{N - 1)(1 - e)(^i + \V2\/S) 

This completes the proof of Lemma 7.8. D 

7.5. Proof of Proposition 7.2. From (P5), one can choose C^ sufficiently 
small such that /" is a density. By (7.1) and (7.6), 



Jt 



x'{g\g-) = Ci5-^^^' --- ,7 ^" d^v 



{{i)^ *w){v )Y 



l2 



-^^^ Jv 9>l,V2/6) ^*" 

Let 

Ju = {v:0<vi< 5, \V2\/S < 6}, Ji2 = {v-0<vi< 6, \v2\/S > 6}, 

J21 = {v:vi> 6, \v2\/S < 6}, J22 = {v.vi > 6, \v2\/S > S}, 

and for i,j = 1,2, 

, _ f {Iv'^u,vAv2-U2)dWs{u)}\ 

By (7.6), Lemmas 7.5, 7.7 and the Plancherel formula (3.9), we obtain 
111 < Ce^^^ I 1 1 ^u,vAy2-U2)dW5{u)\ d,v 

Lemma 7.1 and (P4) imply 



(7.7) / ^u,vA'^2-U2)dWs{u)<Cr]5{v). 

Jv 

Let ci = {N — 1)(1 — ,^) — h. It follows from (7.7), Lemmas 7.7 and 7.8 that 
/i2 = 0((5e"(^i"^)^), /2i = 0(e^('^^'"')^), /22 = 0(<5e-(2-i-i)^). 

Now letting e = min(37r — 26, {N — 1)(1 — ,^) — 26 — 1) > and combining 
the above bounds, we obtain 

X'(5°, 9") < Cld-^'Pihi + I12 + hi + /22) < C5~^^^\-''. 
Choosing 6 = e/logn, we have the desired result. 
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